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Abstract. The following problem, arising from medical imaging, is addressed: Suppose 
^ I that r is a known tetrahedron in R"^ with centroid at the origin. Also known is the 

orthogonal projection U of the vertices of the image (j)T of T under an unknown rotation 
\ (f) about the origin. Under what circumstances can (f) be determined from T and [/? 
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1. Introduction 



The perspective-n-point problem, often abbreviated PnP, is the problem of determining 
the position of a camera from the perspective images of n given points. The problem 
has been widely investigated during the last few decades, using several traditional camera 
models, such as projective (see, for example, [9]), orthographic (see, for example, [H]), or 
|> ■ weak perspective (i.e., scaled orthographic, see [21 [10]), and focusing on various aspects 

O . (such as small values of n). 

While the solution of a specific instance of PnP is often an application of elementary 
I geometry, understanding the configuration space — for example, classifying which config- 

urations admit a given number of solutions — involves challenging nonlinear aspects (cf. 
[6], |T3] and the references therein). Indeed, it was not until recently that Faugere et 
al [6] (partially) classified the configurations for the perspective-3-point problem via the 
discriminant variety, using extensive computations. 

Our point of departure is a paper by Robinson, Hemler and Webber [H], who, motivated 
by an application in imaging, studied the perspective-4-point problem for the orthographic 
cd \ camera model. The problem is as follows. A given tetrahedron T in with vertices 

p^^\ . . . has been transformed by an unknown (direct) rigid motion (p. Also given is 
the image U = {u^^\ . . . , u^^^} of the set of vertices of (pT under a parallel projection onto 
the x|/-plane in an unknown direction w G S*^. The problem is to find (f) and w. 

In [13] it is observed that one may as well take the parallel projection to be the orthog- 
onal projection tt^ onto the xy-plane. It is also noted that then can only be determined 
up to a vertical translation, because such a translation does not change U. Since is the 
composition of a rotation about the origin and a translation, it suffices to determine the 
rotation and the horizontal component of the translation. The authors of [H] make the 
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assumption that it is known which projection comes from which vertex of T, that is, they 
assume that , 4. Under this labehng assumption, they show that 

the rotation and horizontal shift can be determined. 

Our purpose here is to study this problem when the labeling assumption is removed, and 
to provide a systematic foundational study from the viewpoint of nonlinear computational 
geometry (see, for example, [Tl 151 fT^). 

Clearly, the centroid of the vertices of 0T must lie on the vertical line through the known 
centroid of U. From this, we make two conclusions. Firstly, the horizontal shift can always 
be determined, so we may assume that </> is a rotation about the origin. Secondly, if such 
a rotation (j) can be determined when the centroid of T is at the origin, then it can also be 
determined when the centroid of T is located elsewhere. Thus our problem can be stated 
in the following form. 

Suppose that T is a known tetrahedron in with vertices p^^\ . . . ,p^^^ and centroid at 
the origin. Also known is the orthogonal projection U = {u^^\ . . . , u^^^ onto the xy-plane 
of the vertices of the image 0T of T under an unknown rotation (p about the origin. Under 
what circumstances can we determine (j) from T and f/? 

Obviously, if T has nontrivial automorphisms — for example, if T is regular — then (j) 
cannot be uniquely determined. Now let T be an arbitrary tetrahedron in with vertices 
p^^\ . . . ,p^^\ Suppose that the images {{p^^^ +p^^^)/2) and (p [{p^^^ + p^^^)/2j under (p 
of the midpoints of the opposite edges [p^^\p^'^'^] and [p^^^p^'^)] are contained in the ^-axis. 
Then a rotation ip of (pT by tt about the 2;-axis results in a tetrahedron ip(pT whose vertices 
also project onto U. In this case U forms the vertices of a parallelogram in the xy-p\ane, 
so U has a symmetry (rotation by ir about its center). 

These preliminary remarks show that in general (p cannot be determined if T or ?7 
has extra symmetries. A general goal is to understand if it can be uniquely determined 
otherwise, and if not, to find those T and U that do allow (p to be determined. 

The relation between our problem and the one considered in [13] can be made clearer 
if we regard the labels of the vertices of T as having been permuted by an unknown 
permutation a of {1,2,3,4}, so that u^^'^ is the projection of (pp^°'^^'>\ i = 1, ... ,4. Then 
the problem in [13] corresponds to the case when a is the identity. 

In this paper, we deal with both uniqueness and reconstruction. Our focus is on the 
geometry of the problem, in particular, the configuration space of all tetrahedra leading 
(for a given rotation) to the same set of projection points as the original tetrahedron. By 
decomposing this space into the union of the spaces corresponding to the various types of 
permutations involved, we can treat the configuration questions from a linear algebra point 
of view. Then, using some nonlinear symbolic methods, we precisely classify situations 
where the dimension of the configuration space deviates from the expected dimension. As 
a consequence, we are able to prove in Theorem 18.11 that for almost all tetrahedra T in 
with centroid at the origin, there does not exist a rotation other than the identity 
such that 7iz(pT = tt^T. However, the various lemmas that we prove along the way provide 
much more detailed information. 
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Problems such as the one addressed here, involving the retrieval of information about a 
geometric object from data concerning its projections onto lines or planes (or intersections 
with lines or planes), fall under the umbrella of geometric tomography [8]. 

The paper is structured as follows. After the preliminary Section |2l the case when 
the permutation a is the identity is considered in Section |3] from a linear algebra and 
symbolic viewpoint. Then, in Sections IH El El and O we deal with the other cases. 
Finally, in Section |H] we state the main conclusions for our study. 

2. Notation and preliminaries 

As usual, S^~^ denotes the unit sphere and o the origin in Euclidean n-space M". The 
Euclidean norm is denoted by || ■ ||. Unless specified otherwise, Xi will signify the ith 
coordinate of a point x = {xi, . . . , Xn) in M". The unit ball in will be denoted by 5". 
We write [x, y] for the line segment with endpoints x and y. Orthogonal projection onto 
the xy-plane in is denoted by vr^. Given m G S*^, we denote the line through the origin 
parallel to u by /„. The dimension dim A of a set A in M" is the dimension of its affine 
hull. The symmetric group on {1, 2, . . . , n} is denoted by S„. 

Let S0(3) denote the group of rotations about the origin in M.^. An element of S0(3), 
henceforth simply called a rotation, can be specified in terms of a rotation axis (a line 
through o) and a rotation angle. We shall also use the following characterization using 
quaternions (see, for example, [71 Sec. 8.2]). For a quaternion q = a + bi + cj + dk, where 
a, b,c,d & M, the rotation matrix R{q) associated with q is 

. / a^ + b'^ -c^ -d^ 2bc - 2ad 2bd + 2ac 

(2.1) R{q) = ——l 2bc + 2ad - b"^ + - d'^ 2cd - 2ab 

™ \ 2bd - 2ac 2cd + 2ab -b"^ -c^ + d^ 

where ||g|| = \/a? + b'^ + + d^. Conversely, the quaternion q corresponding to a rotation 
with axis in the direction w = {wi, W2, w^) G S*^ and rotation angle a is 

(2.2) q = cos(a/2) + wi sin(a;/2)i + W2 sin(a/2)j + sin(a/2)/c. 

Since a rotation around an axis of rotation lu by angle —a is the same as a rotation around 
/_„ by angle a, we may without loss of generality restrict a to the interval [0, tt]. 

It will be convenient to regard a tetrahedron in M.^ simply as a set of four points in R^. 
Either these points are in general position, in which case they form the set of vertices of 
a full- dimensional tetrahedron in the usual sense of the term, or they are contained in a 
plane and hence lower dimensional. 

Throughout, we consider a known tetrahedron T = {p^^\ . . . , p^^^} in R'^ with centroid at 
the origin. The projection U = {u^^\ . . . , u^^^ of (j)T onto the xy-plane, where (j) G SO (3) 
is unknown, is also given. Then there is an unknown permutation a G S4 such that 

(2.3) 7r,#(^) = 
for i = 1, . . . , 4. 

The case dealt with in [I^, corresponding to a = id, the identity permutation, can be 
viewed as that of a labeled tetrahedron; the projections of the vertices retain the labels, so 
that it is known which point in U corresponds to which vertex of T. 
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If G SO (3) and a G S4, we denote by Tcr{<j)) the family of (possibly lower-dimensional) 
tetrahedra T = {p^^\ . . . ,p^'^^} such that 

(2.4) {<I>P% = pf^'^\ 

for i = 1, . . . , 4 and j = 1,2. When there are two different rotations of T giving rise to 
the same set U of projections onto the xy-plane, we may for our purposes assume that 
one rotation is the identity, and then by (12. 3p . (12.41) holds for some 7^ id. Thus it 
suffices to study the system (12. 4p in order to understand uniqueness issues, and we will 
be interested in the dimension of 7^(0) in various situations. Since each tetrahedron is 
a set of four points in M^, we could also regard a tetrahedron as a point in (R^)^ ~ M}^, 
and thus consider Ta{4>) as a set in M}"^. However, we are assuming that T has centroid 
at the origin, so that 

4 4 

(2.5) J]p«=o = 5^0p«. 

i=l i=l 

Using (12.51) . we may identify T with any three of its points, say the first three, and then 
the equation in (12.41) corresponding to z = 4 is redundant. We shall therefore identify 
7^(0) with the corresponding set in MP, which, in view of (12. 4p and (12. 5p . is actually a 
subspace of M^. Of course each tetrahedron gives rise to not one but 24 points in MP 
(depending on which three of its vertices are selected and in which order), but since we 
are only interested in the dimension of 7^(0), this loss of bijectivity is unimportant. 

Clearly there are only five essentially different cases to consider. There is the labeled 
case when a = id, and if a 7^ id G 5*4, then a is a two-cycle, a direct product of two 
two-cycles, a three-cycle, or a four-cycle. Corresponding to the four latter cases, we can, 
without loss of generality, consider in turn (i) cr = (2,1,3,4), (ii) a = (2,1,4,3), (iii) 
a = (2, 3, 1, 4), and (iv) a = (2, 3, 4, 1). 

3. The labeled case: cr = id 

In the terminology introduced in the previous section, Robinson, Hemler, and Weber 
[H] proved the following result. 

Proposition 3.1. Let T be a full- dimensional labeled tetrahedron in M^ . Then there do 
not exist two different rotations such that the resulting (labeled) projections of the rotated 
vertices ofT onto the xy-plane coincide. Thus the rotation is uniquely determined by the 
(labeled) projection. 

Proof. Suppose two different rotations as in the statement of the proposition exist. Clearly 
we may assume that one is the identity id and denote the other by 7^ id. If the resulting 
projections coincide, then from (12. 4p with a = id, we obtain 

(3.1) (0pW)^.=p«, 
where i = 1, . . . ,4 and j = 1,2. Let 

(3.2) Hj = {x eM^ : (0x)j - Xj = 0}, 
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for j = 1,2. Then p^*-* G -ffi fl if2 for i = 1, . . . , 4. If either of the subspaces Hi and H2 are 
proper subsets of M^, we are done, since T (Z Hi (1 H2 contradicts the assumption that T 
is full dimensional. Otherwise, we have {4>x)j — Xj = for all a; G M'^ and j = 1,2. But 
then (j) fixes the xy-plane and hence = id. □ 

The authors of [H] gave a different proof of the previous proposition, deriving it from a 
reconstruction procedure. For the convenience of the reader, we provide a different recon- 
struction method that can be obtained from the reconstruction result for the perspective- 
3-point problem under weak perspective (see [21 [ID])- Extending the reconstruction algo- 
rithm in [To] to the four-point case works as follows. 

Let be the unknown rotation. We first construct the unique circle C containing the 
known points p'-*-' for ? = 1, 2, 3. We aim to construct the projection E = tt^^C of (pC, an 
ellipse in the xy-plane whose semi-major axis has length equal to the radius of C. The 
known points u^^^ for i = 1, 2, 3 lie on E. For i = 1,2, 3, denote by m^*-' the midpoint of the 
edge of the triangle p^^\ p^'^\ p^^^ opposite to p^*^ and by t*^*^ the other intersection of the 
line through p^*) and m*^*) with the circle C. The corresponding midpoints vr^^m^'^ of the 
edges of the triangle u^^^ , m^^-* , u^^^ opposite to m*^*^ can of course be constructed since this 
triangle is known. Then, for i = 1, 2, 3, the point vr^^t*^*^ can be constructed by elementary 
geometry, since 

_ 77iW|| ||^(») _ 7i^(f)m^''^\\ 

for z = 1, 2, 3. Since vr^^t*^*) lies on E for i = 1, 2, 3, we have constructed six points on E. 
But any five points determine an ellipse, so we can construct E itself. Now E determines 
the circle (pC, up to reflection in the xy-plane and vertical translation, and hence the 
points (j)p^^\ (f)p'^'^\ and are similarly determined. Since T is known, the position of 
(pp'^^'^ is also known relative to 0C, up to a reflection in the plane containing (pC. If T is 
full dimensional, det (p is determined by the points (l)p^'^\ i = 1, . . . , 4, and since det 0=1, 
no reflection is possible. Now we can use the fact that because the centroid of T is at the 
origin, the centroid of 0T is also. This allows 0T and hence (since we are in the labeled 
case) to be completely determined, if T is full dimensional, and up to a reflection in the 
a;|/-plane, if T is contained in a plane. Note that if T is contained in a plane, a reflection 
of T in the xy-plane is of the form ipT for some rotation ip about the origin, so cannot 
be fully determined in this case. 

Recall that we regard the family 7^(0) as a set in and that we are considering the 
case a = id. 

Lemma 3.2. Let (p ^ id be a rotation. Then dim7id(0) = 3, unless the axis of rotation 
is horizontal, when dim7Id(0) = 6. 

Proof. Let T be a tetrahedron with vertices i = 1, ... ,4 and centroid at the origin. 
Identifying T with p^^\ i = 1,2,3 and using (12. 4p with a = id and (12. Sp . we see that 
T G 7Id(0) if and only if 



(3.3) 



(0p»)^._p»=O, 
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for i = 1,2,3 and j = 1,2, a system of six equations in nine variables. Let M be 
the corresponding 6x9 coefficient matrix, where the variables are ordered p^^\ pg^**, 
p"i \ . . . , ^3'^^ and where for i = 1, 2, 3, rows 2z — 1 and 2i of M correspond to the equations 
with index j = 1 and 2, respectively. Then dimTTdl^) equals the dimension of the null 
space of M. 

Since M obviously has rank at most six, we obtain dimTTdl*/*) > 9 — 6 = 3 directly from 
the Rank Theorem. 

Let 
(3.4) 

+ &2 _ c2 - £ 2bc -2ad 2hd + 2ac\ r M 

2hc + 2ad a^-b'^ + c^- d'^ 2cd - 2ab I 1 



The rotation can be represented by the matrix (12.11) with + b'^ + + d"^ = 1, and 
using the latter equation and ( 13. 3p . we can rewrite M as a block matrix, 

A- 1 
M= I A-I 

A-I 



where 

A-I 



— d^ be — ad bd + ac 
be + ad -b"^ — d"^ cd — ab 

Suppose that dim7rd(0) > 3. Then the rank of M is less than six, so all the 6x6 
minors of M vanish. The 6x6 minor corresponding to columns 1, 2, 4, 5, 7, and 8 of M 
is 

— d"^ be — ad 
be + ad -b^ - d^ 

Hence d = 0, which in view of (12. 2p implies that the axis of rotation is horizontal. From 
the geometry it is clear that without loss of generality, we may suppose that this axis is 
parallel to (1, 0, 0), so that b = sin(a/2) ^ and c = 0. But then 



-b^ -ab 



8 det ( T , , V i I = 8d'ia' + b^ + c^ + d^ = 8d^. 



A 



in which case all 4 x 4 minors vanish but not all 3 x 3 minors do so. Then the rank of M 
is three, so dim7rd(0) =9 — 3 = 6. □ 

The geometry corresponding to the previous lemma is as follows. We know from Propo- 
sition [XT] that each member of 7Id(0) is degenerate, and hence contained in a plane. If 
is a rotation about a line not contained in the xy-p\ane, then the only solutions to (13.31) 
are those for which each point pW, i = 1,2,3, is contained in the axis of rotation. For each 
p^^^ there is one degree of freedom, and hence the set of solutions is three dimensional. 
Suppose, on the other hand, that is a rotation by angle a about a horizontal line. Then 
the points p^'^\ i = 1,2, 3, must lie in one of the two planes containing this line and at an 
angle a/2 to the xy-plane. For each p^*^ there are two degrees of freedom, so the set of 
solutions is six dimensional. 
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Note that in the previous discussion about dim7rd(0), the position of p^^'^ is determined 
by the centroid condition fl2.5p . once the positions oi p^^\ i = 1,2, 3, are known. We shall 
use this fact frequently in the sequel without special mention. 

The geometric statements in Proposition 13.11 and Lemma [3^ yield an algebraic corollary. 
To formulate this, let / be the ideal in the real polynomial ring 

R = R[p^^\pS\pi'\p?,...,pt) 

generated by the linear polynomials 

(3.5) 

for i = 1, . . . , 4 and j = 1,2. An ideal generated by linear forms is also called a linear 
ideal. 

Corollary 3.3. For any rotation 7^ id, a positive power of the polynomial 



(3.6) det 



is contained in the linear ideal I. 



/I 1 1 1 \ 

Pi^'^ pf^ pf'' Pi^^ 

(1) (2) (3) (4) 

P2 P2 P2 P2 

, (1) (2) (3) (4) , 

\P3 P3 P3 P3 J 



Proof. By Proposition 13.11 whenever a sequence of points {p^^\ . . . ,p^^^) is a zero of the 
polynomials (13. 5p then the points p^^\ . . . ,p^^'' are affinely dependent, that is, the determi- 
nant (13.61) vanishes. This determinant can be seen as a polynomial in p^^^ 's. Thus, by the 
weak form of Hilbert's NuUstellensatz (see, for example, ^ Section 4.1]), this determinant 
polynomial is contained in the radical ideal rad(/) = {r ^ R : & I for some n G N} . 

We remark that though the NuUstellensatz is a statement over the complex numbers, 
standard Grobner basis theory implies that the determinant is also contained in the real 
linear ideal. Namely, since the ideal is generated by real polynomials, the standard al- 
gorithms for computing a Grobner basis of rad(J) (see, for example, [21 Theorem 8.99]) 
always keep coefficients within the reals and thus provide a real basis for rad(/). Similarly, 
the algorithm for reducing a real polynomial with respect to a Grobner basis generated by 
real polynomials keeps coefficients within the reals. Since for any polynomial in the ideal 
this reduction algorithm yields a representation in terms of the generators, our remark 
follows. □ 



4. One TWO-CYCLE: cr= (2,1,3,4) 

Lemma 4.1. Let (p ^ id be a rotation by angle a and let a = (2, 1, 3, 4). Then dim 7^(0) = 
3, unless (i) the axis of rotation is neither horizontal nor vertical and a = tt, when 
dim 7^(0) = 4, or (ii) either the axis of rotation is vertical and a = n or the axis of 
rotation is horizontal and < a < n , when dim 7^(0) = 5, or (iii) the axis of rotation is 
horizontal and a = tt, when dim 7^(0) = 6. 
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Proof. Let T be a tetrahedron with vertices i = 1, . . . , 4 and centroid at the origin. 
Identifying T with p^^\ i = 1,2,3 and using (12.41) and (12.51) . we see that T E Ta{(l)) if and 
only if 

(4.1) (0p(i)), - pf = 0, (0p(2))^. _ p(^) = 0, and (^p^^^), - pf = 0, 

for i = 1, 2, 3 and j = 1, 2. Let Mi be the 6x9 coefficient matrix of this system of six 
equations in nine variables, where the variables are ordered p^i \ p'^\ p^\ p^i \ . . . ,p^^\ 
and where for i = 1,2, 3, rows 2i — l and 2i of Mi correspond to the equations with index 
j = 1 and 2, respectively. Then dim 7^(0) equals the dimension of the null space of Mi. 
From (14.11) . we have 

( A -I 
Mi= [ -I A 

\ A-I 

where A and I are given by (13.41) . 

Since Mi obviously has rank at most six, we obtain dim 7^(0) > 9 — 6 = 3 directly from 
the Rank Theorem. 

Let Ji be the ideal generated by all 6 x 6-minors of Mi together with the polynomial 
T = a? + h'^ + + (P' — 1. A Grobner basis Gi of Ji with respect to the lexicographic 
ordering a h c d is given by 

Gi = {a^d\ a^cd^ a'^d\ aV, a%cd\ a'd^h^ - - d''), ad\d^ - 1), acd^, a^d\ ahd\ 
ahcd^, alP'd^, r}. 

(This can be found with a variety of standard software. Experts may well prefer a 
different choice, but with Mathematica, it can be done by defining the matrix Mi, 
using Minors [Ml, 6] to generate the 6x6 minors of Mi, adjoining the polynomial 
+ 6^ + + (i^ — 1 to this list, and then using GroebnerBasis [{list} , {a,b, c ,d}] .) 
From this Grobner basis we see that if dim 7^(0) > 3, then a = or d = 0. One can check 
that the rank of Mi is 5, 4, or 3, when a = and 7^ 0, ±1, or when either a = and 
(i = ±1 or a 7^ and c? = 0, or when a = d = 0, respectively. This yields dim 7^(0) for 
cases (i), (ii), and (iii) in the statement of the lemma. □ 

In order to describe the geometry behind the previous lemma, for i = 1,2, let Hi be 
the plane containing j)'-*-' which is orthogonal to the axis of rotation /„, let Ci be the circle 
in Hi containing p^*^ and with center on lu, and let k be the vertical line through p^'^\ See 
Figure [H The rotation takes p'-^^ on li around the circle Ci to the point on I2 and 
also takes p^'^^ on I2 around the circle C2 to the point cpp^'^'^ on li. The angle of rotation is 
of course the same in each case, and we also have 

since the planes Hi and H2 are parallel and so intersect h and I2 in equidistant pairs of 
points. It follows that Ci and C2 have equal radii and hence p^^^ and p^^^ are the same 
distance from /„. 

If lu is neither vertical nor horizontal, then vr^Ci and 11^02 are ellipses with their centers 
on Tizlu- If TT^Ci 7^ TizC2, these two ellipses intersect in two points, namely vr^p^^^ = vr^^p^^^ 
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Figure 1. Geometry behind Lemma 4.1. 



and TT^p^^^ = TXz(l)p^^\ which are reflections of each other in tXzIu- Moreover, the angle a 
of rotation must be strictly between and vr. Since Ci and C2 must intersect the vertical 
lines li and I2 through these two points, there is only one degree of freedom in choosing 
the position of each of p^^^ and p^'^\ The point p^^^ must lie on /„, allowing a further degree 
of freedom, so there are a total of three degrees of freedom, as in the first statement of 
the lemma. If vr^Ci = 'n:zC2i then Ci = C2 and a = n. In this case there are three degrees 
of freedom in choosing p^^\ after which the position of p'-^-* is determined, and one in 
choosing p^^\ This corresponds to case (i) in the statement of the lemma. 

Suppose that lu is the 2;-axis. Then Ci and C2 are possibly different horizontal circles 
and a = TV. There are three degrees of freedom in choosing p^^^ and one each for p^^^ and 
p^^\ since the latter point must lie on the 2;-axis. This situation is included in case (ii) in 
the statement of the lemma. 

Finally, suppose that lu is contained in the xy-plane. If < a < vr, then Ci = C2 is a 
vertical circle. There are three degrees of freedom choosing p^^\ after which the position 
of p^^^ is determined, and two degrees of freedom in choosing p^^^ (which must lie in the 
plane containing lu and having angle a/2 with the xy-plane). Again, this is included in 
case (ii) in the statement of the theorem. If a = tt, then Ci and C2 are possibly different 
circles lying in the same vertical plane. There are three degrees of freedom in choosing 
p*^^) and then one degree of freedom for p^^^ , since it can lie anywhere on the vertical line 
through (pp^^\ The point p'-^^ lies in the vertical plane containing /„, allowing two degrees 
of freedom in its choice. Thus there are six degrees of freedom in total, as in case (iii) in 
the statement of the lemma. 



5. Two TWO-CYCLES: o-=(2, 1,4, 3) 

Lemma 5.1. Let (p ^ id be a rotation by angle a and let a = (2, 1, 4, 3). Then dim 7^(0) = 
3, unless (i) the axis of rotation is horizontal and < a < n, when dim 7^(0) = 4, or (ii) 
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the axis of rotation is neither horizontal nor vertical and a = tt, when dim 7^(0) = 5, or 
(iii) the axis of rotation is horizontal and a = ir, when dim 7^(0) = 6, or (iv) the axis of 
rotation is vertical and a = tt, when dim 7^(0) = 7. 

Proof. As in the proof of Lemma [4.H T G 7^(0) if and only if 

(5.1) (0p(^)),- - pf = 0, (0p(2)) . _ pW = 0, and {(f)p^^^)j + pf^ + pf + pf = 0, 

for z = 1, 2, 3 and j = 1,2, where in the third equation we have used fl2.5p to write p^^^ 

in terms of p^j \ i = 1,2,3. Let M2 be the 6x9 coefficient matrix of this system of six 
equations in nine variables, under the same convention as in the proof of Lemma I4.1[ 
From (15.11) . we have 

I A -I 
M2 = I -I A 

\ / / A + I 

where A and / are given by (13.41) . 

The matrix M2 has rank at most six, so dim 7^(0) > 9 — 6 = 3. 

Let J2 be the ideal generated by all 6 x 6 minors of M2 together with the polynomial 
r = + 6^ + + — 1. A Grobner basis G2 of J2 with respect to the lexicographic 
ordering a>~b>~c>~d is 

G2 = {a^d\d^ - 1), a^cd{d^ - 1), a^c^d, a%d{d^ - 1), a%cd, a\b^ - c^)d, a^cd, a%d, t}. 

From this we see that if dim 7^(0) > 3, then a = or d = 0. (Note that d = ±1 implies 
a = b = c = 0.) One can check that the rank of M2 is 5, 4, 3 or 2, when a 7^ and d = 0, 
or when a = and d 0, ±1, or when a = d = 0, or when a = and d = ±1, respectively. 
This yields dim 7^(0) for cases (i), (ii), (iii), and (iv) in the statement of the lemma. □ 

The geometry behind the previous lemma is straightforward using the analysis given 
after Lemma [4.11 and bearing in mind the centroid condition (12.51) . We omit the details. 
Note that case (iv), when is a rotation by tt about the z-axis, was already mentioned 
in the introduction. 

6. Three-CYCLE: a = (2,3,1,4) 

Lemma 6.1. Let (p ^ id be a rotation by angle a and let a = (2, 3, 1, 4). Then dim 7^(0) = 
3, unless (i) the axis of rotation is horizontal, in which case dim 7^(0) = 4, or (ii) the 
axis of rotation is vertical and a = 27r/3, in which case dim 7^(0) = 5. 

Proof. As in the proof of Lemma [4. H T G Ta{<f)) if and only if 

(6.1) (0P^^^)i - pf = 0, (0p(^))j - pf = 0, and (0p(^^)j- - pf = 0, 

for i = 1, 2, 3 and j = 1, 2. Let M3 be the 6x9 coefficient matrix of this system of six 
equations in nine variables, under the same convention as in the proof of Lemma 14. 1[ 
From (16.11) . we have 

/ A -I \ 
M3 = A -I \, 
\-I A J 
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where A and / are given by (13.41) . 

Then M3 has rank at most six, so dim 7^(0) > 9 — 6 = 3. 

Let J3 be the ideal generated by all 6 x 6 minors of M3 together with the polynomial 
T = a"^ + b"^ + + — 1. A Grobner basis G3 of J3 with respect to the lexicographic 
ordering a>~b>~c>~d is 

G3 = {c/2(4^2 _ 3^2^ ^^^4^2 _ ^^(4^2 _ 3^^ ^(^2 ^ ^2^^ 

It follows that if dim 7^(0) > 3, then either d = or d = ±a/3/2 and b = c = (and 
hence a = ±1/2). One can check that the rank of M3 is then either 5 or 4, respectively. 
This yields dim 7^ for cases (i) and (ii) in the statement of the lemma. □ 

Again, we comment on the geometry behind the previous lemma. If the axis of rotation 
lu is horizontal, there are three degrees of freedom in choosing p^^K Then p*^^) and p'-'^^ 
must lie in the vertical plane H containing p'-^^ and orthogonal to lu- Moreover, p^^^ must 
lie in the line obtained by rotating the vertical line through p^^^ by —a around (so 
that (pp^^^ and p^^^ have the same projection on the xy-plane). This is another degree of 
freedom. Similarly, p^^^ must lie in the line obtained by rotating the vertical line through 
p(3) _Q, around lu (so that (pp^"^^ and p^^^ have the same projections on the xy-plane). 
But p*-^-* must also lie in the vertical line through (f)p^^^ (so that ^p*-^-* and p*-^-* have the 
same projections on the xy-plane). This means that p*^^-* is determined by the positions 
of p*-^-* and p^^^ and so there are only four degrees of freedom in this case. 

If the axis of rotation is the 2;-axis and the angle of rotation a = 27r/3, there are three 
degrees of freedom in choosing p^^^ and a further one degree of freedom for each of p^"^^ 
and p^^\ since their heights may be different from that of p^^-* and only their horizontal 
positions are determined. Thus there are five degrees of freedom in all, corresponding to 
case (ii) in the statement of the lemma. 

In the general case, there are three degrees of freedom in choosing p^^\ after which the 
other points are determined. Indeed, p*-^-* must lie on the vertical line /i, say, through 
(j)p^^\ and p^^^ must lie on the vertical plane through the line obtained by rotating li by 
a around lu (so that 0p(^) and p^^^ have the same projection on the xy-plane). Moreover, 
p(3) must also lie on the line obtained by rotating the vertical line through p^^^ by —a 
around lu (so that (pp^^^^ and p^^^ have the same projection on the xy-plane). Hence, in 
the general case, p^^-* is determined and consequently also p^'^\ 

7. FoUR-CYCLE: (7= (2,3,4,1) 

Lemma 7.1. Let (p ^ id be a rotation by angle a and let a = (2, 3, 4, 1). Then dim 7^(0) = 
3, unless (i) the axis of rotation is neither horizontal nor vertical and a = it, in which 
case dim 7^(0) = 4, or (ii) the axis of rotation is vertical and either a = tt /2 or a = tt, 
in which case dim 7^(0) = 5. 

Proof. As in the proof of Lemma [4. 1^ T G 7^(0) if and only if 

(7.1) (0p(^)),- - pf = 0, (0p(2)),- - pf = 0, and (0p(3))j + pf^ + pf + pf = 0, 
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for i = 1, 2, 3 and j = 1,2, where in the third equation we have used fl2.5l) to write Pj 

in terms of p^j \ i = 1,2,3. Let M4 be the 6x9 coefficient matrix of this system of six 
equations in nine variables, under the same convention as in the proof of Lemma I4.1[ 
From fl7.ll) . we have 

I A -I 
M4= A -/ 
V / / A + I 

where A and / are given by (13.41) . 

Since M4 has rank at most six, we obtain dim7^(</)) > 9 — 6 = 3. 

Let J4 be the ideal generated by all 6 x 6 minors of M3 together with the polynomial 
T = a? + h'^ + + d? — 1. A Grobner basis G4 of J4 with respect to the lexicographic 
ordering a y b >~ c >~ d turns out to be 

= {a\2d^ - 1)2, a'c(2rf2 - 1), a%{2d^ - 1), a\b^ + c^), a{d^ - l){2d^ - l)^ 
ac{2d'^ - 1), ab{2d^ - 1), a{b^ + c^), r}. 

Consequently, if dim 7^(0) > 3, then either a = or d = ±l/v^ and b = c = (and 
hence a = ±1/^2). It can be verified that the rank of M4 is 5 if a = and d ^ 0,±1 or 
4 if either a = and d = ±1 or a = ci = ±1/a/2 and b = c = 0. This yields dim 7^ for 
cases (i) and (ii) in the statement of the lemma. □ 

Regarding the previous lemma, suppose that the axis of rotation is not horizontal or 
vertical. If the angle of rotation a = n, there are two degrees of freedom for choosing p*^^^ 
in the vertical plane containing ly^, after which p^"^^ and p^^^ can be chosen anywhere in 
the vertical line containing p^^\ making four degrees of freedom in all. This corresponds 
to case (i) in the statement of the lemma. 

If lu is the ^-axis and a = 7r/2 or a = n, there are three degrees of freedom for choosing 
p^^\ After this the horizontal positions of p'-^^ and p^^^ are determined but their heights 
are arbitrary, giving a total of five degrees of freedom. This deals with case (ii) in the 
statement of the lemma. 

It remains to explain the generic case. It is clear that there are no solutions when 
lu is the z-axis unless a = 7r/2 or a = n, and it is easy to see that if lu is horizontal, 
then we have three degrees of freedom in choosing p^^\ after which the other points are 
determined. Suppose, then, that is neither vertical nor horizontal and < a < tt. For 
i = 2,3,4, let n^*^ = {p^^^ +p^^^)/2 be the midpoint of the edge [p^^\p^^^]. Using (12. 5p . we 
obtain 

= (p(l) +p(2) _p(3) _p(4))/4^ 

(7.2) n(3) = (p(i) - p(2) + p(3) _ p(4) ) 

= (p(l) _ p(2) _ p(3) + p(4))/4. 

Notice from (17.21) that the points n^*^ determine T via the equations p^^^ = rS^'^ +n^^'^ + rS'^\ 

p(2) = „(2) _ „(3) _ ^(4)^ p(3) _ _^(2) ^ ^(3) _ ^(4)^ ^(4) _ _^(2) _ ^(3) ^ ^(4)_ ^Yisll 

therefore focus on the degrees of freedom in specifying n^*-*, z = 2, 3, 4. 
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To this end, suppose that 

(7.3) #«=p(*+i)+/i,e3, 

for 2 = 1, . . . , 4, where /ij G M, indices are taken modulo 4, and 63 denotes the unit vector 
in the direction of the positive 2-axis. From fl7.2p and fl7.3l) . we get 

0n(2) = -n(^) + (/ii + /i2)e3/2, 

(7.4) 0^(3) = _^(3) + (^^+^^)g^/2^ 

<^^(4) = ^(2) + + /i4)e3/2. 

For z = 2, 3, 4, let Ci be the circle contained in a plane orthogonal to n, with center c'-*-' 
on and containing n*^*). We claim that the position of c'-'^^ on lu alone determines that of 
rS^\ and hence there is only one degree of freedom in choosing n^^\ To see this, note that 
— C3 is the circle orthogonal to u, with center — c^^^ on and containing —rS^\ Let u-^ 
be the plane through the origin orthogonal to u and let vr denote the parallel projection 
in the direction 63 onto M"*". Then txCj, and 7r(— C3) are circles in u-^ with centers ttc^^^ 
and 7r(— c*^^)) on 7r/„. See Figure [2l By the second equation in (17.41) . 7r(0n^^^) = 7r(— n^'^^). 
Since 7r(0?T,^^^) lies on 7rC3, we see that 7rC3 and 7r(— C3) intersect at 7™^^-* and 7r(— n^'^-'). 
Now vrn*-^^ and n^—rS'^^) lie on a line through the origin, so this line is orthogonal to vr/^. 




Figure 2. Geometry in the plane u after parallel projection vr. 

Therefore we have vrn'-^'' = A(m x 63)/!! u x 63 1| , for some real A, because tto*^^-* is orthogonal 
to both u and 63. The angle between vm^^^ and T:{—rS^^) = ir{(j)n^^^) at ttc^^^ is a, so 

A = llvm^^^ll = IIttc^^^II tan(a/2). 

It follows that once the position of c^'^^ on is specified, we know nc^'^^ and therefore A 

1 /2 

and hence 7rn^^\ We also know the radius (||7rc*-^^|p + Hvrn*^'^-' |p) of vrCs, which equals 
that of C3. From this and Tm^^\ the position of n''^-* is determined. This proves the claim. 
Next, we consider n^^-* and n^^\ The first and third equations in fl7.4p tell us that 

(7.5) 7r(0n(2)) = 7r(-r2(^)) and 7r(0n(^)) = vm^^). 

Identify M"*" with the complex plane C in such a way that irl^ is the real axis. Then since 
TTc*^^^ and vrc^^^ lie on 7r/„, they are real. Let u = exp{—ai). Then by (17. 5p . we have 

(7m(2) - 7rc(2))a; = -nn^^^ - ttc^^) , 
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If cj^ + 1 7^ 0, we can solve the linear system (]7.6p for vrn'-^^ and vrn''^^ in terms of vrc^^-* and 
7rc^^\ Therefore once the positions of c*-^^ and c*-^-* on are specified (for which there are 
two degrees of freedom), we know vrc'-^^ and 7ic^^\ hence irn^'^^ and 7in^^\ As above, this 
allows the radii of the circles C2 and C4 to be determined, and then n^"^^ and n*-^-* are also 
determined. Finally, if + 1 = 0, then a = ±7r/2. Then it is easy to see that for (17. 5p 
to hold, we must have 7rC2 = nC^. This means that is vertical, which is not the case. 




Figure 3. A four-cycle example (perspective view). Dotted lines through 
the common projections i = 1, ... ,4, of pairs of vertices of the two 
tetrahedra onto the xy-pla.ne are vertical. 



Example 7.2. For a specific example of the four-cycle situation, let u = (l/-\/2, 0, 1/V2) 
and let be the rotation around by 7r/3, so that has matrix 

3/4 -v^S/S 1/4 
V^STS 1/2 • 
1/4 3/4 

Let = (-2, -3 + 76, 16 - 3^6), p^^) = (1^ 3 _ 4^5^ _i9 + 3^5), = (^2, -3 + 
3\/6, 8 — 3-\/6), and p(^) = (-1,3,-5 + 3 a/6). Then it is easy to check that the tetrahedron 
r = {pW,...,p(4)}isfull dimensional, and = (1, 3-4^6, 13-3^6), ^ (2,-3 + 
sVQ), -20 + 3^6), = 3^ II _ 3^5)^ g^j^^ 0p(4) ^ ^_2^ .3 + _4 + 3^), The 
projections onto the xy-pl&ne give the set U = {u^^\ . . . , u^^^, where u^^^ = (—2, — 3+-\/6), 
^(2) = (1,3 _ 4^6), = (2,-3 + 3^6), and u^^'> = (-1,3)}. See Figure 1 for an 
illustration of this example. 



DETERMINING A ROTATION OF A TETRAHEDRON FROM A PROJECTION 



15 



8. Main results 

Theorem 8.1. For almost all tetrahedra T in with centroid at the origin, there does 
not exist a 7^ id G S0(3) such that 7rz(pT = tt^T. Indeed, the exceptional set constitutes 
a finite union of subspaces, each of dimension at most seven, in M^. 

Proof. Let T be a tetrahedron in M.^ with centroid at the origin, and suppose that 7^ 
id G SO (3) is such that vr^^T = vr^T. Then there is a (Jq G S4 such that (12. 4p holds with 
a = ao and hence 

TGr.„(0)c U TM- 
(J6S4 

By our conventions and Lemmas 13. 2[ 14. I5.H I6.H and I7.H the latter set is a finite union 
of subspaces of M^, each of which has dimension at most seven. Therefore this set is of 
zero Lebesgue 9-dimensional measure and the theorem is proved. □ 

Example 8.2. We claim that a specific example of a full- dimensional tetrahedron satis- 
fying Theorem 18.11 is T = {p^^\ . . . where 



P 



(1) 



;i,0,0), p^'^ = (1,1,0), p(3) = (2,1,2), and p(^) = (4, -2, -2). 



To see this, observe that 

(8.1) lb^'^f = l, lb^'^f = 2, |b(')|r = 9, and Ib^f = 24, 
while 

(8.2) ||vr.p«ir = l, ||7r.p(^)|r = 2, \\nj'^\' = 5, and ||7r.p(^)f = 20, 
By (12^ . we have 

for i = 1, . . . ,4. Comparing (18. ip and (18. 2p . we see that the only possibility is that a = id. 
Since T is full dimensional, our claim follows from Proposition 13.11 

From a practical point of view, perhaps the most important observation is that there 
are only 24 ways to label the points in a tetrahedron T in M^. If T is full dimensional, 
then, for any particular such labeling, can be reconstructed by the method of Section [31 
or that of Robinson, Hemler, and Webber [T3], or symbolically using standard software, 
yielding at most 24 solutions for the rotation 0. 

We close with a remark illustrating how the uniqueness issues are refiected within sym- 
bolic reconstruction methods. While for each fixed permutation a, there is a unique solu- 
tion for reconstructing a full-dimensional tetrahedron (since fixing the permutation allows 
Proposition 13 . 1 1 to be applied), there may be more than one solution for lower-dimensional 
tetrahedra. For example, consider the tetrahedron T with vertices p^^^ = (—1, 0, 1), p^^-* = 
(0, 0, 0), = (0, 0, -2), and p(^) = (1, 0, 1), with n^T = U = {(-1, 0), (0, 0), (1, 0)}. Let 
a = (2, 3) be the one-cycle that interchanges 2 and 3. Solving symbolically, we obtain 
the Grobner basis {d,c,b^ — 6, a6, + 6^ — 1}, which yields four distinct solutions for 
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{a,b,c,d). However, these only result in two rotation matrices, 

1 \ / 1 \ 

10 and 0-10. 
001/ \00-l/ 

Of course, these correspond to the two possible rotations (p such that tt^^T = U, namely, 
the identity and the rotation by n about the x-axis. 
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